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We study the propagation of non-relativistic bound states moving at constant velocity across a 
homogeneous thermal bath and we develop the effective field theory which is relevant in various 
dynamical regimes. We consider values of the velocity of the bound state ranging from moderate to 
highly relativistic and temperatures at all relevant scales smaller than the mass of the particles that 
form the bound state. In particular, we consider two distinct temperature regimes, corresponding 
to temperatures smaller or higher than the typical momentum transfer in the bound state. For 
temperatures smaller or of the order of the typical momentum transfer, we restrict our analysis to the 
simplest system, a hydrogen-like atom. We build the effective theory for this system first considering 
moderate values of the velocity and then the relativistic case. For large values of the velocity of 
the bound state, the separation of scales is such that the corresponding effective theory resembles 
the soft coUinear effective theory (SCET). For temperatures larger than the typical momentum 
transfer we also consider muonic hydrogen propagating in a plasma which contains photons and 
massless electrons and positrons, so that the system resembles very much a heavy quarkonium 
' in a thermal medium of deconfined quarks and gluons. We study the behavior of the real and 

^ |. imaginary part of the static two-body potential, for various velocities of the bound state, in the hard 

thermal loop approximation. We find that Landau damping ceases to be the relevant mechanism for 
dissociation from a certain "critical" velocity on in favor of screening. Our results are relevant for 
understanding how the properties of heavy quarkonia states produced in the initial fusion of partons 
in the relativistic collision of heavy ions are affected by the presence of an equilibrated quark-gluon 
plasma. 

(N 
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^ , 

CN , I. INTRODUCTION 

' When matter is immersed in a thermal medium many of its properties change. In principle, no strictly stationary 
^-^ ' bound state exists, because interactions with the particles of the medium lead to a finite lifetime for all states 
(including the ground state). This is equivalent to a broadening of the energy levels, i.e. an imaginary part of the 
energy eigenvalues, which depends on the density and on the temperature of the medium. 

Of particular interest is the case in which the bound state moves with respect to the thermal medium. The first 
k> ' experimental investigations and theoretical developments of this system were done in condensed matter physics 
; I From the analysis of atoms moving across a plasma, it has been shown that a number of phenomena may take place. 
First of all the Debye screening of the Coulomb potential depends on the relative velocity between the bound state and 
the plasma. Moreover, the propagation of a bound state through the medium produces a fluctuation of the induced 
potential which leads to a density trail. Finally, the moving particle loses energy and is eventually stopped by the 
plasma. 

A renewed interest in the properties of bound states moving in a thermal medium arose in recent years due 
to the advent of high energy heavy-ion colliders. In particular, one is interested in understanding whether some 
modifications of the properties of heavy quarkonia (HQ) states produced in the early stage of the heavy-ion collision 
can be a signature of the presence of a deconfined plasma of quarks and gluons. In their pioneering work 2], Matsui 
and Satz showed that the Debye screening of the color interaction between two static heavy quarks may lead to the 
dissociation of HQ in a thermal medium. This effect should be experimentally detectable by the suppression of the 
corresponding yields. The suppression of HQ states means that the yield of HQ observed in heavy-ion collisions is 
smaller than the yield of HQ one would obtain multiplying HQ production rates in p-p collisions by the number 
of nucleons participating in the collision and taking into account the normal nuclear absorption; see e.g. [3] for a 
brief review. The first study of moving HQ was then performed in [4], in which the dependence of the Debye mass 
on the velocity of propagation of the heavy quarks with respect to the quark-gluon plasma (QGP) was determined. 
Subsequent analyses have confirmed the effect and studied the formation of wakes in the QGP |5l-[lO|. 
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One may wonder whether the drift of bound states is important, as it is the case for heavy flavors. Indeed, 
measurements of heavy flavor production in PHENIX [llj via single electron measurements results in a large V2, 
which suggest that there is significant damping of heavy quarks while they travel across the fireball. Therefore, in 
heavy-ion collisions the thermal bath expansion may drift the heavy quarks in a phenomenon similar to advection in 
normal fluids. This pict ure has also received support from microscopic calculations of heavy quark diffusion in the 
quark-gluon plasma |12l |. However, we expect that the drag of a heavy quarkonium is less important than that of a 
heavy quark. This is because an isolated heavy quark has a net color charge while a heavy quarkonium at distances 
larger than its radius is colorless. Hence, in general, we expect that the HQ states produced in the early times of the 
collision will not be comoving with the thermal medium, and, therefore, our calculations will be relevant for them. On 
the other hand HQ states produced through recombination are expected to roughly comove with the thermal bath. 
This is because both heavy quarks have been drifted by the QGP before recombining in a HQ. 

Suppression of the J/^* was first observed at the CERN SPS However, in contrast with the naive Debye 

screening scenario, further experimental investigation of the J/^E" yields at PHENIX led to the observation of a 
strong suppression at forward rapidity rather than at mid-rapidity. Recently, there have been efforts in studying this 
problem with the use of non-relativistic effective field theories (EFTs) [ijl - llq . The EFT techniques are very useful 
for problems that have different energy scales, as is the case of HQ in a thermal medium. Using these techniques it 
has been shown that, at least in perturbation theory, the dissociation of bound states is due to the appearance of an 
imaginary part in the potential [H, [l^ [13] . 

In the present paper we study how a moving thermal bath affects the properties of bound states. One of the 
points is to assess whether the results obtained in a static medium are modified when considering the relative motion 
between the bound state and the thermal medium. We consider the simplest systems, hydrogen-like atoms moving at 
a constant velocity, v, across a homogeneous thermal medium. We study two different cases, the first one corresponds 
to temperatures smaller or of the order of the typical momentum transfer, the second one corresponds to temperatures 
larger than the typical momentum transfer. We always assume that the temperature is much smaller than the mass 
of the particles forming the bound state. 

In the first case we restrict ourselves to the hydrogen atom. We consider separately temperatures smaller than the 
typical momentum transfer and temperatures of the order of the typical momentum transfer. In both temperature 
ranges we provide the matching procedure and evaluate the energy shifts and decay widths for the stationary states of 
the system. We build the effective theory for this system first considering moderate values of the velocity and then the 
relativistic case. We show that for large values of the velocity, a new separation of scale occurs and a different EFT 
must be constructed for dealing with bound states. In this case the separation of scale is such that the corresponding 
EFT resembles some aspects of the soft coUinear effective theory (SCET) p^ . 

In the second case, namely, for temperatures larger than the typical momentum transfer, we also consider muonic 
hydrogen. Since the mass of the particles forming the bound state is much larger than the mass of the particles in the 
thermal bath, this system resembles very much a heavy quarkonium in a thermal medium of deconfined quarks and 
gluons. As a consequence, this part of the present work is of direct relevance for understanding how the properties of 
heavy quarkonia states produced in the initial fusion of partons in the relativistic collision of heavy ions are affected 
by the presence of an equilibrated quark-gluon plasma. We study the behavior of the real and imaginary part of the 
two-body potential for various values of the velocity of the bound state with respect to the thermal bath employing the 
hard thermal loop (HTL) approximation. Regarding the real part of the potential we reproduce known results, and 
extend them to higher speeds. The imaginary part of the potential is calculated for the first time. We demonstrate 
that screening overtakes Landau damping as the dominant mechanism for dissociation at a certain critical velocity. 

This paper is organized as follows. In Section |ll] we introduce some general remarks about the propagation of 
particles in a thermal bath. In Section [IIII and Section ITVl we study the hydrogen atom moving at moderate velocities 
and ultrarelativistic velocities with respect to the medium respectively. In Section |V] we study the real and imaginary 
part of the static potential of muonic hydrogen in a moving thermal bath for temperatures larger than the typical 
momentum transfer. The results of this last section are directly applicable to the HQ case. Finally, we present our 
conclusions in Section IVll 

II. GENERAL FRAMEWORK 

In our study we shall employ a reference frame in which the bound state is at rest and the thermal medium 
moves with a velocity v. We choose this frame because it facilitates the application of non-relativistic effective field 
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theories, in particular, Non-Relativistic QED fi^ (NRQED) and potential NRQED ^ (pNRQED)i. These EFTs are 
extremely convenient to handle the three different scales of non-relativistic systems at vanishing temperature [2l|, , 
and have already proved useful to analyze these systems in a static thermal bath, in which additional scales occur 
[23} . They allow one to organize the calculations in such a way that only one scale is taken into account at each 
step, which, together with the use of dimensional regularization, makes computations much easier. 

We shall assume that the plasma (or black-body radiation) is in thermal equilibrium at a temperature T. Since we 
are considering the reference frame in which the plasma is moving with a velocity v the particle distribution functions 
are given by 

where the plus (minus) sign refers to fermions (bosons). In the reference frame where the thermal bath is at rest 
pt^kf^^i^, while in a frame where the plasma moves with a velocity v we have that 



/3^ = ^(l,v) = ^, (2) 



where 7 = l/Vl — v'^, v = |v|, is the Lorentz factor. This frame has been successfully used in the past, for example 
in 1241 ■ Studying a bound state in a moving thermal bath is akin to study a bound state in non-equilibrium field the- 
ory [25| ; in that case the Bose-Einstein or Fermi-Dirac distribution functions are substituted by a general distribution, 
which in our case will be the boosted Bose-Einstein or Fermi-Dirac distribution functions reported in Eq. ([Ij. 

The vector v brings in the problem a number of complications. First of all, it breaks rotational invariance. Second, 
when V gets close to 1 new scales are induced, which has serious implications for the use of EFTs. For instance, if 
we are working with an EFT in which we have integrated out all scales larger than /z, we can no longer argue that if 
/i ^ T the Lagrangian of this EFT is not affected by the temperature. This is because the Boltzmann suppression 
is not only controlled by T, but rather is a non-trivial function of T and v. In order to illustrate this point, let us 
analyze the distribution functions in Eq. ([T]) in more detail. 

We begin with a thermal bath consisting of massless particles. Taking into account that in non-equilibrium field 
theory the collective behavior always enters through on-shell particles or antiparticles, we have (in the case of particles) 
that 

where k = |k| and 9 is the angle between k and v. The distribution functions in Eq. ([T} can now be written as 

where we have defined the effective temperature 

T,^{e.v) = -^ -. (5) 

1 — V cos u 

Intuitively, the dependence of the effective temperature on v and 9 can be understood as a Dopplcr effect. Indeed 
Eq. ([5]) is analogous to the change in the frequency of light caused by the relative motion of the source and the observer. 
Consider a particle in a thermal bath of radiation moving with velocity v. From the point of view of the particle 
it will see in the forward direction blueshifted radiation and in the backward direction redshifted radiation. This 
corresponds, respectively, to an effective temperature which is higher in the forward direction than in the backward 
direction. Therefore, the effective temperature corresponds to the temperature of the radiation as measured by the 
moving observer and by a minor abuse of language we shall talk about blueshifted and redshifted temperatures. Notice 
that analogously to the relativistic Doppler effect, the effective temperature in the transverse direction, i.e. in the 
direction corresponding io 9 — 7r/2, is redshifted. 

For w ^ 1, one has that Tce{9,v) ~ T for any value of 9 and one single scale T controls the Boltzmann factor in 
Eq. (HI. However, for v close to 1, the values of Tcs{9, v) strongly depend on 6*, which gives rise to an interesting case 



Although the Lagrangian of NRQED is known for an arbitrary reference frame, most of the developments have been carried out in the 
rest frame of the bound states. 
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for an EFT analysis. In order to proceed further, it is convenient to use light-cone coordinates. We choose v in the z 
direction and define 

fc+ = /cq + ^3 and fc_ = fco — fca . (6) 

Then, we have that 
where 



Therefore, in light cone coordinates, it becomes explicit that the distribution function actually depends on two scales, 
T-|_ and T_ . For any value of v it is clear that Tj^ > T > r_ and moreover T4. correspond to the highest temperature 
measurable by the observer, while r_ corresponds to the lowest temperature measurable by the observer. 

For small values of w, T+ ~ T_ and the shift in the temperature in the forward and backward directions are negligible. 
In this case no further separation of scale is needed and one has a single temperature scale, T, as mentioned before. For 
u ~ 1, we have T_|_ 3> T^, namely, two well separated temperature scales, which must be properly taken into account 
in our EFTs. Note that configurations with light-cone momenta such that ^ or /e_ ^ T_ are exponentially 
suppressed. Then we can separate the remaining configurations in two regions (in light-cone momenta): 

• A collinear region, corresponding to fc+ ^ and fc_ < r_ . 

• An ultrasoft region, corresponding to <C and k_ < T_ . 

The existence of these two regions has to be taken into account in the matching procedure between different EFTs. 
In this paper we shall analyze two different situations in which v is close to 1, the case nie ^ ^ 1/r ^ T_ 3> E 
and the case r+ ~ me > 1/r :^T_:^ E. 

We would like to remark that although in Eq. ([3|) we have assumed for simplicity that the particles of the thermal 
bath are massless, our approach applies to a plasma that consists of both massless and massive particles. Note that 
our discussion, from Eq. ^ on, holds independently of what the dispersion relation of the particles in the thermal 
bath is. If some particles in the plasma have mass M ^ T, we know that they are exponentially suppressed in the 
thermal bath, and this must be true in any reference frame. We can easily verify it by substituting k_ = (A/^+k^)/fc+ 
in Eq. ([7]), and by noticing that its minimum is attained at /J'^fc^ = M/T, which confirms that for M ^ T thermal 
effects due to these particle can indeed be neglected. 



III. HYDROGEN ATOM AT MODERATE VELOCITIES 



In the present section we shall assume that the velocity is moderate, say v < 0.5, so that ~ T ~ T_, and 
separately study the cases T <C 1/r and T ^ 1/r (r is the size of the bound state, and hence 1/r of the order of the 
typical momentum transfer). For simplicity, we also assume that the proton is infinitely heavy. 

As we have explained in the previous section, in a thermal bath at a temperature T ^ M, particles with mass M 
are exponentially suppressed independently of the value of v. In particular, if M ~ me indicates the mass of electrons 
and positrons in the plasma, then these particles are irrelevant in our analysis. In this range of temperature the hard 
thermal loop effects will not appear and the doubling of degrees of freedom plays no role [ij] . This implies that only 
the transverse photons are sensitive to thermal effects and, hence, the leading order interaction, namely, the Coulomb 
potential, will not be modified. 



A. The T < 1/r case 

As a starting point we consider the pNRQED Lagrangian at vanishing temperature (we use the form given in Eq. 
(6) of dil). We will be able to evaluate the corrections to the binding energy En and to the decay width r„ due to 
the thermal bath up to the order ruea^- There are two different diagrams that contribute at this order. The first one 
is the tad-pole diagram which is given by 

V ■/ ~~^eJ (27r)^-i cI^'-'^mI ~ 1 " 3^' ^' 
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where the soHd lines represent the atom propagator and the wavy Une corresponds to the photon propagator. In the 
integral we set the number of dimensions D — A, because the integral is convergent and fc'' corresponds to the loop 
momentum. Notice that the contribution of this diagram is independent of v because the loop integral has no indices 
and no external momentum enters in it. This is in fact true for any tad-pole diagram of this kind. Therefore one can 
read the result from the v = case [ij] . 
Next we consider the "rainbow" diagram, 



= 7T lim 



(10) 



is the momentum operator of the electron, \r) symbolizes an eigenstate of the Coulomb Hamiltonian of energy Ej- , 
and Iij{q) is defined as follows: 



d^k sm - k^) 



(27r)'D-i el/5^'=f'l -Iq-ka+ie 



~fc2~ 



Since the only two independent tensors are 5ij and ViVj, one can use the decomposition 



where 



and then 



^(4 



2 



and 



„2 



Upon substituting Eq. (jlll) in the expressions above, we find that 

d^k 6iki-k^) i ^ ^ f d^k 5{kl-k^) 



A ^ 



and 



B 



(v-k)2 



(11) 
(12) 
(13) 

(14) 
(15) 



The computations of these integrals is done in Appendix [A] The results for the imaginary and real part of A are, 
respectively, given by 



q 

47r2 



log 



(—) 

\ kl J 



2v 



log 



1-v 



ttTVT 



\q\v 



-3 log 




and 



(16) 



Sttv 



los 



The imaginary part of B turns out to be given by 



^B{q) = 



127r2 



lot 



while for the real part of B we have 



tVT 



2v^ 



1 -e 



dXX"^^^ 



\2ttT^1 - 



dt 1 



1 



(17) 



(18) 



(19) 



The expressions above can be computed numerically for any value of the parameters. The thermal corrections to the 
energy and the decay width (for arbitrary angular momentum) are given by the following expressions: 



Irn 



3m 



— lim > 



;n|p>')3/,,(po-i?.)(rb^|n). 



(20) 



6 



and 

<5r„z™ = ^ lim y^{n\p'\r)m,,{po- Er){r\p^n) , (21) 

which in general will depend on the relative velocity v. Analytical expressions for T ^ £' and for T E, where 
E ~ En, the binding energy scale, are derived below. 



1. TheT^ E case 



For T ^ E the leading contribution to the integrals in Eqs. (|17l) and Eqs. (|19p can be analytically 

determined and upon substituting the corresponding expressions in lij we find that 



8ttv 



1-v 



-OiE), 



(22) 



for the real part of lij and 



47r2 



^ log 



2ttT 



„ 1 / 27rT 



11 1 

^+3^2^ 



1-t) 



(23) 

for the imaginary part of lij . 

It is interesting to observe that, in this limit, the terms non-local in q (Bethe-log type) coincide with those obtained 
for vanishing velocity. Hence, all dependence on v is encoded in an anisotropic potential and kinetic term. This result 
will also serve as a cross-check of the 1/r ~ T computation that will be carried out in the next section. 

In order to obtain the energy shift and the decay width from Eqs. (1^^ and P5)) we consider separately the S-wave 
states and states with non-vanishing angular momentum. We display below slightly more general results, which hold 
for an ion of charge Z as well. 

In the S-wave states the expected value of any tensor operator {n\Oij\n) oc Sij, therefore terms proportional to P[j, 
can be ignored (because P[j is traceless) and we find that 



SEn 



2a 



AZa^ |0„(O)p 



2v ^\l-v 



+ 1-7 



J2\{n\p\r}\^{En-Er)\og 



2t:T 



where (/)„(0) is the wave function at the origin. The corresponding change in the width turns out to be given by 



(24) 



2Z^a^TVl - , /I 
<jr„ = — log 



3n u 



1 



(25) 



States with non- vanishing angular momentum are more difficult to deal with. It is convenient to decompose lij by 
the tensors Sij and instead of Pfj and P[j . For the real and imaginary parts of lij we find respectively 



16ttv 



(26) 



and 



- — 



7+g 



1 3 

4^2 8v 



^ 1 



2^^(") 



(27) 



where 



(28) 
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To determine the energy shifts and the decay widths we fix v in the z-direction, and make use of the following identities 



(v.p)2 



Air 



Yi„ 



and 



{n\[[H,^r■ p], v p]|n) 





1 









n)(rz™|F20|i1m), 



(29) 



(30) 



where Yim are the spherical harmonics. Note that \n) is used as a short-hand notation for \nlm), where n is the 
principal quantum number, I the orbital angular momentum and m its z-component. With these expressions we 
obtain the general forms of the shifts of the energy levels 



5E. 



2a 



nlm 



3me 37rmg ^ 
Z^a^{2m\l0){2l0m\lm 



2nmlalnH{l + + 1) 
and the corresponding widths are given by 

'1 + v 



Y,\{n\p\r)\'{E^-Er)log 



-El 



\En — Er 



Z^a^TVT 



2 log 



1-v 



1 



log 



1 + v 
1 - V 



{2l00\l0){2l0m\lm) , 



(31) 



(32) 



where {lml'm'\l"m") are the Glebsch-Gordan coefficients (normalization and sign conventions are as in [2f 
It is interesting to observe that the decay widths and ([5^ decrease as the velocity increases. 



2. TheT <^E case 

For temperatures T E the coefficients A and B simplify and upon replacing their expressions in Eq. (jl2p we 
obtain that 

^«<') = (T^ + ^'^«) + K?) ■ '''' 

When this expression is used in the evaluation of the energy shift in Eq. and of the decay width in Eq. (|2ip , we 
need to calculate 

r 

for X — s,p. Note that only the term proportional to Pfj contributes because P[j is traceless. Thus, the contribution 
from the rainbow diagram in the limit T <^ E for states with vanishing angular momentum is given by 



3m„ V E^ 



and is independent of the velocity v. Hence, the dominant contribution of the rainbow diagram above cancels the 
contribution of the tad-pole diagram, which we have seen to be independent of v as well. Then the thermal corrections 
in this case are very suppressed, of the order 0{a^T'^ / E^), like in the case of the thermal bath at rest. 



B. The r~l/r case 



In the temperature regime T ^ 1/r, the temperature is high enough so that its effects must be taken into account 
already in the matching between NRQED and pNRQED, namely it affects the potential. Therefore, several diagrams 
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are modified by the presence of the temperature. However, like in the calculations for the thermal bath at rest, there 
are only four diagrams that give a relevant contribution. All other diagrams give contributions that either vanish or 
can be shown to cancel out by local field redefinitions. We schematically analyze the relevant diagrams below. Then, 
as a cross check, we compare the calculations for T ~ 1/r in the limit of low temperature, with the results that we 
derived in the previous section for E ^ T in the limit of high temperature, and we find agreement. 

The four diagrams that must be taken into account in the matching procedure between NRQED and pNRQED are 
the following: 



The tad-pole diagram which comes from the term in the NRQED Lagrangian and gives the contribution 

(36) 



2me 

3mp 




This diagram is quite similar to the diagram we evaluated in the previous section, but now the solid line 
corresponds to the electron field instead of the hydrogen atom. As already discussed, this diagram is independent 
of the velocity v. 



The rainbow diagram is given by 



1 1 \ 

vT + 



where the solid line corresponds to the electron field and the wavy line corresponds to the photon. In order to 
have a consistent EFT we have to expand 

-2?^-^, (38) 



and therefore the contribution of the rainbow diagram can be written as 



me 



where 



and 



- J (2^)^-ifc3 ^^^(1--!^) _ ^ (^^^ p ) ' (40) 



jD-li 



(2^)^-'fce^7fc^^'"'^'-l 
For the time being we do not evaluate these integrals; as we shall clarify soon we only need to evaluate Ti 



• The thermal correction of the Coulomb potential corresponds to the diagram 



-^^^ =^^T.,, (42) 



where the tensor Tij is the same as in Eq. (I40[) and p and p' are the momenta of the incoming and outgoing 
electrons, respectively. The solid thick line here corresponds to the ion propagator and the solid thin line is the 
electron propagator. 
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• The last diagram to consider is the relativistic tad-pole, that is the same as the previous tad-pole diagram, but 
now the vertex comes from the term in the NRQED Lagrangian. The contribution of this diagram is given 

by 




6m5 



——■n'n^E- 



(43) 



where the tensor Rij is defined in Eq. (|41]) . Notice that the term on the right hand side of Eq. ([39|) cancels the 
corresponding contribution of this diagram. Therefore, the sum of the rainbow diagram and of the relativistic 
tad-pole diagram is independent of Rij . 

• The remaining non- vanishing diagrams give contributions analogous to the ones of Eqs. (36), (37) and (42) 
of Their net effect can be shown to be zero by local field redefinitions, like in the case of the thermal bath 
at rest. 



Now, let us consider how these diagrams combine. In particular, we would like to obtain the pNRQED Lagrangian 
that matches all these terms. By inspection we find that the pNRQED Lagrangian is given by 



6C 



pNRQED 




2me 



4a7r Za 
ml Ixi - X2I 



3rne ^ ^ 6mg 
T„-W^(t,xi)V^'7^(i,Xi), 



(44) 



which as already noted does not depend on Rij. We evaluate Tij in a similar way as was done for lij in the T ^ E 
case. The result in the M S subtraction scheme is 



T = 



67r2 



log ( ^ ) - log2 - log(27r) + 7 + -1 ( 1 + -i^ ) log 



T 



8v 



3v^ 



1 

4^2 



Pjv) VjVj 

87r2 



(45) 



where we have decomposed Tij in terms of 5ij and ViVj /v^ , and p{v) is defined in Eq. 



Upon making a local field 



redefinition to remove the term with a time derivative in (1441) . one can identify the corrections to the potential, in a 
similar way as it was done in Eq. (45) of The final form of the thermal correction to the pNRQED Lagrangian 
reads 



5L 



pNBQED 



+ 7 + 



— n ^ 

8w V 3w2 



log 



6m|! 
'1 - 



j d^xi (fx2N\t,X2)N{t,X2) 



log 



2ttT 



-log 2 



1 



4w2 



^3(xi -X2) + 



V'^(i,xi)V'(i,xi), (46) 



where Vc{r) is the Coulomb potential at vanishing temperature. With simple modifications this Lagrangian can be 
put in a form so that we have an atom field instead of an electron and a nucleus field (see 14] for more details). 

Now that we have computed the corrections to the pNRQED Lagrangian for the case T ^ 1/r we also need to 
compute the contribution from the ultrasoft scale with this Lagrangian. These contributions can be computed from 
the tad-pole and the rainbow diagrams as in Eqs.® and ([TU]), where the Bose-Einstein distribution function can be 
expanded because we are now in the case T ^ E, and therefore 



(47) 



Upon substituting the expansion above in Eq. © we find that the contribution of the tadpole diagram vanishes in 
dimensional regularization, because it has no scales (it is independent of the external momentum). The rainbow 
diagram gives a contribution similar to the one in Eq. ()10p , with the replacement lij — >■ Jij , where 



^Jij — 



8ttv 



K log 



1 



lot 



pp.. 



(s) 



2v 



(48) 
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and 



In order to be consistent with pS)) . the MS scheme has also been used here to remove the UV divergences. The 
thermal corrections to the energy levels and to the decay width coming from the soft and the ultrasoft scales in the 
case T ^ 1/r are respectively given by 

airT'^ na^T'^ ^Zo? /,/«\ 3/, l\ 1 



- ("1^5 + ^ Jim^E(-b»»^^.(^'o - ^r)(r\^\^) , (50) 



and 



'5r„i^ = ^ lim V(n|p»3fiJy(po-Sr)(r|p^|n). (51) 

Note that the dependence in the correction to the Darwin term in ([50)1 is canceled out by the /i dependency of Jjj (g) 
in Eq. (gHl)- Note also that, upon substituting in (|5T|) . the expression for the decay width reduces to that of ([5^ . 
Furthermore, the thermal corrections to the binding energy above in the limit of low temperature coincide with (p4)) 
and (|3ip . Therefore, the limit of low temperature in the case T \/r agrees with the limit of high temperature in 
the T ^ E case. 



IV. HYDROGEN ATOM AT RELATIVISTIC VELOCITIES 



When the bound state moves at high speed with respect to the thermal bath one has to take into account that 
the effective temperature measured by the bound state in the forward direction is blueshifted and that the effective 
temperature in the backward direction is redshifted. In particular one has that ^ T_, and therefore and r_, 
which we have defined in Eq. (|8]), are two well separated energy scales that must be properly taken into account in 
the EFT. In particular it is possible that and T_ are in two distinct energy ranges. In this case the analysis of 
the system differs considerably with respect to the case of the thermal bath at rest. We shall study two different 
situations of this sort: the first one corresponds to the case ^ 1/r ^ T- ^ E and the second one to the case 
T+ ~ me ^ l/r ^ T_ 3> Recall that we are assuming that T ^ Wg and hence, as we have already stressed in 
Section ini we can neglect positrons and electrons in the thermal bath. 



A. The T+ ~ 1/r > r„ > case 



Since r+ <C Wg we can use NRQED at vanishing temperature as the starting point, but we would like to integrate 
out also the l/r scale in order to construct the pNRQED for this situation. As in the soft coUinear effective theory, it 
is convenient to split the photon field A^^ into two different components, a coUinear one A™' that takes into account 
photons with fc+ ^ and A;_ ^ T_ and a ultrasoft one A^^ that takes into account photons with k+ ^ ^ T_ 
(or smaller). Notice that both types of photons have virtualities A that fulfill (1/f)^ ^ A; this means that neither 
of these two types of photons have to be integrated out in the matching between NRQED and pNRQED. We carry 
out the matching below using an electron field and a nucleus field in pNRQED (rather than an atom field) , as it was 
done in [2l|. 



1. Matching between NRQED and pNRQED for colUnear photons 



In the matching procedure between NRQED and pNRQED we have to determine the effective vertex between non- 
relativistic electrons and coUinear photons. In pNRQED the interaction of non-relativistic electrons with coUinear 
photons cannot be given by the minimal coupling diagram shown in Fig. [T] for kinematical reasons, as we argue next. 
Let us call p the momentum of the incoming electron and k the momentum of the incoming coUinear photon. The 

2 

non-relativistic electron in pNRQED must have po — ^ l/r, because rrie/r is precisely the typical virtuality of 
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FIG. 1: Minimal coupling in pNRQED between an electron and a collinear photon. The solid lines represent the electron 
propagators, the zigzag line represents the collinear photon. This diagram is not kinematically allowed because when a non- 
relativistic electron interacts with an ultrasoft photon its virtuality changes by a quantity of the order nie/r. Then, the outgoing 
electron cannot be described within pNRQED. 



FIG. 2: Leading-order term for the interaction of electrons with collinear photons in pNRQED. 



the electrons that have been integrated out in the matching between NRQED and pNRQED. However, if the photon 
is collinear, namely ^ 1/r, then the virtuality of the outgoing electron is of order me/r, in contradiction with the 
fact that electrons with such a virtuality do not appear in pNRQED. 

Therefore, the interaction between electrons and collinear photons in pNRQED is given at leading order by 4-point 
processes as the one presented in Fig. [S] The matching procedure is outlined in the following equation 




(52) 

The NRQED diagrams on the left hand side have to match the pNRQED diagram on the right hand side. This part 
of the pNRQED Lagrangian involving collinear photons is given at the required order in Appendix IC II The part of 
the pNRQED Lagrangian involving ultrasoft photons only is the same as in the case with the thermal bath at rest. 



2. Computation using pNRQED 

Since we have determined the pNRQED Lagrangian, we can now calculate the contribution of thermal collinear 
photons to the self-energy of the hydrogen atom [from the part of the Lagrangian reported in Eq. (jCl[) ]. We shall use 
the Coulomb gauge and therefore only spatial components contribute. Moreover, the condition V • A = for collinear 
photons means that A^{x) ^ A^{x) and we only need the first and second terms of Eq. (|C1[) . The contribution of 
collinear photons to the self-energy of the hydrogen atom is given by 

^ + - (53) 




where the zig-zag line corresponds to a collinear photon and the solid line represents the hydrogen atom. As already 
noticed in the discussion of Eq. the tad-pole diagram is not sensitive to the relative motion between the bound 
state and the thermal bath, because no external momenta enter into the loop. The corresponding shift to the energy 
levels is given by 

3me 

As we shall see soon, the contribution coming from thermal collinear photons is the dominant one, however it does 
not depend on the quantum numbers of the state and therefore it cannot be seen in the emission spectra. 

There are two one-loop contributions of ultrasoft photons to the self-energy of the hydrogen atom. The tad-pole 
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contribution is given by 



ie"^ f d^k S{k+k^ - kl) _ ie^ 17j_2 f°° dk^ f°° 



. 

and we find that in dimensional regularization this integral vanishes. The second contribution is due to the rainbow 
diagram 

2 

/^^^ lim V(n|/|r)^..bo-i?.)(rb^'|n), (56) 



d^k 6{k+k^ kj) f hkj 



where 

- [ 

(27r)^-l^i|i^ + ^|_^ V ) q-\{k+ + k^ 

There are three different integration regions that contribute to this integral 

• the region with k^,k^ ^ r_ , 

• the region with k^,k^ ^ 9 , and 

• the region with k^ ^ q and fc_ ^ q{T^/T-^) . 

Note that k+/T+ <C 1 in all the regions. Evaluating the contributions of each region (see Appendix IB]) and putting 
them together we find that 

K,, = aP^^ + &PP. , (58) 

where P^^ and Pf^ are defined in Eq. and 



^« = f^/^log(^^ (59) 



47rVl + i' \1 — t; 



and 



Note that in the v ^ 1 limit, the coefficients a and 6 are equal to the coefficients A and B reported, respectively, in 
Eqs.([Tni), (H?!) and Eqs.dH]), (HH) in the limit q < T_. Therefore, in the same limit, we have that Kij -> Uj. The 
thermal corrections to the energy and decay widths due to the ultrasoft photons can be written as 

5KL = 4 lim V("IP»S^''J-(Po-i?r)(r|p^"|n), (63) 



and 



^^nim = K lim ^(n|pV)3?i^,,(po - i?.)(rb^"|n) . (64) 



For S-wave states we obtain the following expressions for the energy shifts: 

1 fl-v\\ 2a ^ 

2^°HTT^jj^^"^^'^"'P''" ^"^^^"n 2.T 



4Zq2 / 1 /l-w\\ |(/)„(0)|2 2a /|S„-Sr|\ 
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and decay widths 



"^^"O" 3n2 



1 



■log 



1 - V 



(66) 



For states with non- vanishing angular momentum / we find that 



se: 



'jilni 



2a 



Y,\{n\p\r}\'{E^-Er)\og 



~Ei 



lEn — Er 



Z^a^{2m\l0){2l0m\lm) 
67rm2a3Z(/ + i)(; + l) 



(67) 



and 



nlm — 



3n2 



1 + v 



loe 



l + v\ / fl + v 
- log 



1-v 



1-v 



-3 {2m\l0){2l0m\lm) 



(68) 



The total thermal width is given by the ultrasoft contribution reported in Eq. for S-wave states, or in Eq. 
for states with non-vanishing angular momentum. In order to obtain the total thermal energy shift the collinear 
contribution given in Eq. (|54l) must be added to the ultrasoft contributions given in Eq. (|65p for S-wave states, or in 
Eq. (|67p for states with non- vanishing angular momentum. Note that the latter turns out to be totally independent 
of the velocity. 

Note that the decay widths (|66)) and (|68|) are decreasing functions of the velocity, like in the moderate velocity case. 
Furthermore, the results above agree with those of Sec. IIII Al in the w — > 1 limit. 



B. The r+ ~ > 1/r > r_ > £ case 



We shall now consider a highly relativistic hydrogen atom immersed in a thermal bath at a temperature T ^ 1 jr. We 
shall assume that the relative velocity between the hydrogen atom and the thermal bath is such that the temperature 
in the forward direction is blueshifted to the electron mass, that is Tj^ ^ mg, while in the backward direction the 
effective temperature is redshifted to 1/r ^ r„ ^ E. The effective temperatures and T_ are now very well 
separated scales, therefore this situation is specially suitable for the use of EFT. 

In the construction of the effective theory we start with QED at vanishing temperature, because rrie 3> T. However, 
the existence of collinear photons must be taken into account in the matching between QED and NRQED. On this 
aspect the matching procedure is akin to the one in SCET. We shall schematically describe this matching procedure 
below. Regarding collinear photons, they have a virtuality of order (l/r)^, and they must be integrated out when 
matching from NRQED to pNRQED. Finally, the contributions of ultrasoft photons are calculated in pNRQED. In 
this case pNRQED does not include collinear photons, which already have been integrated out. The interaction with 
ultrasoft photons is exactly the same as in the previous case. Their contribution is given by exactly the same diagram 
as in Eq. (15^ and therefore one obtains the energy shifts reported in Eq. (1551) . for S-wave states, and in Eq. ([S71) . for 
states with non- vanishing angular momentum, and the widths are given by the same expressions reported in Eq. (j66p , 
for S-wave states, and in Eq. ([68)) . for states with non- vanishing angular momentum. 



1. Matching between QED and NRQED for collinear photons 

In QED, when a non- relativistic electron absorbs a collinear photon, it turns into a relativistic electron. This means 
that the NRQED Lagrangian cannot have this kind of 3-body interaction (a similar argument was used in Section 
IIV A 1[) . Hence, in NRQED the interaction with non- relativistic electrons has to be a 4-body interaction. 

In this case, there is the additional complication that on the QED side we have bispinors, while in NRQED we have 
only spinors. This can be solved using the non-relativistic projector. The matching equation takes the form 




(69) 

where Z is the wave function renormalization of NRQED that depends quadratically of the momentum. The result 
for this matching is given in Appendix IC 21 
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To match NRQED with pNRQED we have to integrate out colhnear photons. The contribution of colhnear photons 
to the self-energy is given by 




3me \ 2™^ J ' 



(70) 



where we have used the Lagrangian of Eq. (|C6I) in the Coulomb gauge. Note that in this gauge thermal effects are 
only due to the spatial components of and A3 <C A±, and we only need the terms proportional to ci, C2, cn, and 
C13 reported in Appendix IC 21 This diagram was already evaluated in the case of the thermal bath at rest. Since 
tad-pole diagrams are unaffected by the motion of the thermal bath, the result remains the same, and the only effect 
is a constant energy shift in the pNRQED Lagrangian, which amounts to the following shift of the effective mass of 
the electron: 

Snie = . 71 

Regarding heavy quarks, the net effect of colhnear gluons is a very tiny shift of the heavy quark mass by an amount 
of SniQ ^ asT'^/mq, which is irrelevant for the stability analysis of heavy quarkonia. Much more important for the 
stability analysis is how the Coulomb potential changes at high temperatures and this will be studied in the following 
section. 

V. THE STATIC POTENTIAL OF MUONIC HYDROGEN IN THE RANGE T > 1/r 

In muonic hydrogen the proton is orbited by a muon and the bound state consists of two heavy particles. Since the 
muon is about 207 times heavier than the electron, muonic hydrogen is much more compact than standard hydrogen 
and the energy levels of the system have about 207 times the energy of standard hydrogen. Muonic hydrogen 



is investigated in order to have high precision measurements of the proton properties 27[, mainly by Lamb shift 
measurements [2^ . The study of muonic atoms is also important for muon-catalyzed fusion processes [2^ , which are 
under experimental investigation at RIKEN 'so'] and Star Scientific [sT'l. 

The study of muonic hydrogen in a thermal bath with 3> T 3> rrie is akin to the study of HQ states in the 
quark-gluon plasma with mp ^ T ^ Aqcd 3> ruq, where mg is the mass of the heavy quark and rriq {q = u, d, s) 
is the mass of light quarks ^23i] . The reason is that in both cases the temperature is on the one hand much smaller 
than the masses of the particles that form the bound state and on the other hand much larger than the mass of the 
particles in the thermal bath. There are thermally excited electrons and positrons in the QED plasma and thermally 
excited light quarks in the QGP which can modify the Coulomb interaction between the two heavy particles. Actually, 
we shall assume that the temperature is high enough so that we can neglect the masses of the light particles of the 
plasma. 

Apart from the modification of the static Coulomb potential, the propagation of a particle in the medium produces a 
fluctuation of the induced potential which leads to a variation in the density of the plasma. If the plasma behaves as a 
liquid the moving bound state can produce a wake. These effects were first analyzed in condensed matter physics (see 
e.g. [H) and then studied in the context of heavy- ion collisions [3-[Toj and in strongly coupled = 4 supersymmetric 
Yang-Mills plasmas 

In the present section we study the modifications to the leading-order potential between two heavy sources in 
relative motion with respect to the thermal bath at a velocity v. We evaluate the potential in the HTL approximation 
assuming that the temperature of the plasma is T 3> 1/r. The real part of the potential is screened by massless 
particles loops, both in QED and QCD (the only difference between QED and QCD in our results is, apart from 
trivial color factors, the value of the Debye mass tud)- The real part of the potential between a quark and an 
antiquark moving in a thermal bath was first computed in the HTL approximation in [4| and then more recently 
in Q. Recent perturbative calculations [l7| at vanishing velocity have pointed out the importance of the imaginary 
part of the potential. So far its effect has not been taken into account in a moving thermal bath. 

In the Coulomb gauge the potential is obtained by the Fourier transform of the longitudinal photon propagator, 

An{k) = i[A«(fc) + A^(fc) + As(fc)] , (72) 

for fco ^ |k|, where Aij(fc) and Ayi(fc) are respectively the retarded and the advanced propagators and A5'(fc) is the 
symmetric propagator. For a bound state comoving with the thermal bath, it is enough to compute the retarded 
self-energy in the rest frame of the thermal bath and then using 

A*j,{k) = AAik) , (73) 
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and 

As{k) = [l + 2/(|fco|,T)]sgn(fco)[Afl(fc)- A^(fc)], (74) 

one can determine the potential. In the expression above, f{k(),T) is the distribution function of the longitudinal 
photons in the thermal bath. However, the last relation does not hold for a bound state moving through a thermal 
bath j25i] , and must be substituted by the following one: 

where = 7(1, v) is the 4-velocity. Thus, in order to determine the propagator one has to evaluate the self-energies 
IlR{k,u) and Ils{k,u). 

The retarded self-energy Ilji{k,u) was computed in 4] and here we only show the result in the reference frame 
where the bound state is at rest^ 

Un{k,u) = a{z) + -^, (76) 
1 — ti^ 



where z = , ^^"'^^ — ^ Q is the angle between k and v, and 

Y 1 — sin^ 



= f -l)fl„(i±i±^)). ,77) 

b(z) = (a{z) - rai,(l - z') (l - 1 1„ (^i±i±|^^ ) . (78) 

Regarding the symmetric self-energy of the longitudinal photons Ils{k,u), the computation is similar to the one 
done for the retarded self-energy in Q. Consider the full symmetric self-energy tensor 11^^^,. It obeys the Ward identity 

fc''n^,=0, (79) 

and is symmetric, 

n^. = . (80) 

Then, it must have the following structure: 

''•^ = Hi f g^'' - ^) + , (81) 



where Hi and 112 are two scalars and 



is the component of u^^ orthogonal to fc^. Since 11^^^ is a tensor, we can determine the values of Hi and 112 in any 
reference frame, and it is convenient to consider the comoving frame, i.e. the frame in which the thermal bath is at 
rest. It is useful to define the tensor 

p.^ ^ L^^^ - g.^ + ^ ] , (83) 



2 V k^-{k-uYy 

where 

ki^^kf" -{k- u)u>' (84) 



In Q there is a misprint in the first line of Eq. (8), in which the global sign must be the opposite. 
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is the component of orthogonal to w^, and k']_ = k"^ — {k ■ u)'^. It is clear that P^^^u^ = P^^k^ = 0, and therefore 
P'^" projects four- vectors in the direction orthogonal to and fc'^. By means of Eqs. (15^ and (IMl) it is easy to show 
that P^jyU^ = Ppi/fcjl = 0, as well. Then we have that 

p^'^n^, = p''''5p.ni - -Hi , (85) 

and in the comoving frame one has that the only nonvanishing components of P^i^ are 

p.. ^ 1 ( _ ^) , (86) 



2 



and then in this frame 



1 / A- ' ki \ 

P-n^, = 2 (<5- - — j , (87) 

which is precisely the transverse component of the photon (gluon) self-energy in the Coulomb gauge. This quantity 
has been computed for vanishing velocity in |25| and we find that 



^(fc • U)2 - fc2 y (fc.y)2_fc2 

The scalar quantity u^u'^II^j^ has a simple interpretation in the comoving frame, where it turns out to be given by 
ngg. Then we have that 

\{k-u)^-k^ jy' \ {k-uY-k^j ^ ^{k -uY-k-^j 

and we can now compute the symmetric self-energy in the frame where the muonic hydrogen is at rest and the thermal 
bath is moving with a velocity v. In this frame we have that 

n, i27rm2,T(l-w2)3/2(i + i^cos2 6l) 
I - |k|(l - w2 sm 6')'^/2 

We have all the necessary quantities to construct A5 using Eq. (1751) a nd then the propagator in Eq. (17^ . When the 
limit ti — )■ is taken, we obtain for A5 the same result as in refs. [l4 [TEI. [T7|. 

In the previous discussion we have not distinguished between the case of moderate velocities and the case of 
relativistic velocities, as we did in the hydrogen atom calculations. This is motivated by the fact that the results 
found in Section IIVI are identical to the ones that can be deduced by taking the u — )■ 1 limit of the results of Section 
mil However, it is interesting to sketch how the computation of An would be carried out in light-cone coordinates 
for t; ^ 1. One would start with the NRQED Lagrangian that includes the interaction with collinear photons of 
Eq. (jC6l) . Since the virtuality of the collinear photons is of order and T ^ 1/r, they can be integrated out before 
evaluating the potential. At leading order, this gives rise to an energy shift similar to the one reported in Eq. (j71l) . 
In the light sector of the NRQED Lagrangian there are also interactions between soft photons and collinear electrons 
and integrating out collinear electrons one obtains the HTL Lagrangian. If the scale T- is much smaller than 1 /r one 
should consider its effects in the ultrasoft photons. However, these photons can only give subleading contributions by 
means of loop corrections. From now on we consider that the distinction between moderate and relativistic velocities 
is not essential and will not be done. 

From the Fourier transform of the An propagator we have determined the real and imaginary parts of the potential 
in the HTL approximation. The potential is anisotropic, and in Fig. [3] we display the plots of the real (upper panels) 
and imaginary (lower panels) part of the potential for = 0, w = 0.55, and v = 0.99 respectively. We consider two 
directions: the first is along the direction of movement of the thermal bath (right panels) and the second one is along 
the direction orthogonal to the thermal bath (left panels). We plot only positive values of r because the potential is 
symmetric for r — >■ — r. We normalize the real part of the potential to amo, which describes the typical strength of 
the potential in the v — case. The imaginary part of the potential is normalized to aT. With these normalizations 
the displayed shapes hold both for muonic hydrogen and for heavy quarkonium. 

Regarding the real part, we observe that for v ~ 0.55 it is very similar (in fact, the curves overlap to a large extent) 
to the real part in the v = case, being Debye screened at a distance of order mo and not very asymmetric. For 
t) ~ 1, however, although the real part of the potential remains Debye screened at roughly the same distance, it 
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10 
r m 



FIG. 3: Real (upper panels) and imaginary (lower panels) parts of the potential between a quark and an antiquark moving 
with velocities v — 0, 0.55, 0.99 with respect to the thermal bath. Right (left) panels correspond to the direction parallel 
(perpendicular) to the velocity of the thermal medium. 



develops a rather large anisotropy. Indeed, an oscillation is observed in the direction of motion, which leads to the 
formation of a wake in the plasma. 

Concerning the imaginary part of the potential for v ^ 0.55, it is not very asymmetric and remains very similar 
to the V — case. It monotonically increases and keeps the same pattern until v ^ 0.9. From that velocity on the 
imaginary part decreases and the anisotropy grows (see the v = 0.99 curve). In the direction parallel to v one has 
that there is a mild increase with respect to the v = case for r ~ 4m£), and an oscillatory behavior at larger r is also 
displayed. In the direction orthogonal to v there is an enhancement of the potential at short range and a decrease 
at large distances. Note that the imaginary part of the potential vanishes in the origin for any value of v and in any 
direction. 

In Fig. m we plot the contour lines for the real and imaginary parts of the two-body potentials between two 
heavy particles with opposite charges, for various values of the velocity. We focus on the short distance regime (the 
normalization for the real part is slightly different from the one in Fig. [3] because we have inverted the sign in order to 
match the normalization of ^J]), because the distances relevant for dissociation are in the range rmo ^ 1 [l3|- Various 
plots of the contour lines of the real part of the two-body potential were also shown in [3] (see also jg]) and it can be 
seen that we obtained exactly the same results. The contour lines for the imaginary part of the two-body potential 
are reported for the first time here, and we observe that an important anisotropy exists even at short distances. 

Let us next estimate the dissociation temperature in a way similar to refs. [l^ [23j. If we assume that the typical 
momentum transfer k is larger than the velocity dependent screening mass m'jj{v, 9) ^ \I\-ji{k, u)\, we obtain that at 
k ^ 6 3 T-v/l — the real and the imaginary parts of the potential have the same size'^. For moderate velocities no 
qualitative change is expected with respect to the w = case. However, for v close to 1, fc becomes small and it is not 
guaranteed that the screening mass can be neglected in front of it. Indeed, we find that the screening mass remains 



^ This is so for generic 0, meaning d ■k/2. For 9 ~ 71/2, we obtain k ~ e3T/(l — v'^)^^'^. 
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FIG. 4: Contour lines of the two-body potential between a quark and antiquark moving with respect to the thermal bath with 
four different values of the velocity: « = 0, 0.5, 0.9, 0.99. The left panels correspond to the real part of the potential, while the 
right panels correspond to the imaginary part of the potential. The real part of the potential has been normalized to amn, 
while the imaginary part of the potential has been normalized to aT. 
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finite for v close to 1 ^. This means that the typical k for which the real and imaginary parts of the potential have 
the same size is smaller than the screening mass. Since a screened potential only supports bound states of a typical k 
larger than to_d(u, 9), we conclude that at relativistic velocities, unlike the case of moderate velocities, the dissociation 
occurs due to screening (i.e. at the scale Td ~ mec), as originally proposed by Matsui and Satz rather than due 
to Landau damping [ij, [3)[l3]- This can also be qualitatively understood from our plots. For v large and increasing, 
we see from Fig. [3] that the real part of the potential increases whereas the imaginary part decreases. Therefore, 
from some v on, the real part of the potential dominates over the imaginary part and one can find the approximate 
wave functions of the system by solving a standard Schrodinger equation with a real potential. The decay width may 
be then calculated in perturbation theory by sandwiching the imaginary part of the potential between those wave 
functions. The wave functions of bound states go to a vanishing value at the distance where the real part of the 
potential becomes flat. From Fig. 2] it is also clear that the real part of the potential at short distances becomes 
steeper at increasing v. On the one hand this implies that no bound state exists from a certain velocity on. On the 
other hand it implies that when bound states still exist their wave functions are increasingly localized close to the 
origin. Since the imaginary part of the potential goes to zero at the origin, it follows that the decay width of such 
states is also going to zero at increasing v. 

Vc at which screening overtakes Landau damping as the dominant mechanism for dissociation, by equating eT to 
e^/^T{l — v'^y^^ above. We obtain Vc ^ \/l — ae'^/^, where a is a numerical factor of order one ^. A quantitative 
study of all these issues may be carried out by numerically solving the Schrodinger equation with the full (complex) 
potential, for instance along the lines j36t- 38j . This is however beyond the scope of this paper. 

Regarding the real part of the potential, we find it interesting to compare the results reported in Fig. [3] with the 
recent results obtained for super Yang-Mills theory using AdS/CFT [3J|. For this theory it was stated that the 
potential could be approximated with a Yukawa potential and the dependence on the velocity encoded in a screening 
length that depends on v and 6 as follows: 



where h{v, 6) is a function that is almost constant for any v and 0. The expression above does not give a good 
approximation of the potential in the HTL approximation. In particular the Debye screening for v = 0.99 is strongly 
dependent on the angle 0. If we try to fit the exponential behavior of the longitudinal and transverse directions, we 
find a screening length which is about a factor 2 larger in the longitudinal direction with respect to the transverse 
direction. This can also be inferred from the fact that the scale of to|,(i;,6') must be given by |nfl(fc,M)|. We obtain 
in the case d — 1 



Be aware that the 9 above is the angle between the velocity and the momentum transfer, whereas the 9 in (pij) is the 
angle between the velocity and the relative position. In any case, the expression (|92p shows that at ultrarelativistic 
velocities a strong anisotropy for real space potential is expected, as confirmed by our figures and discussed above. 

Regarding the oscillatory part of the potential one might wonder whether it is due to the weak coupling approxi- 
mation. However, as shown in [lo| for the potential produced by a single charge, in the HTL resummation approach 
the oscillations of the potential are larger than in the HTL approximation. Moreover, one would naively expect that 
with increasing coupling the wakes should be larger than in the weak coupling approximation. 



The EFT theory for the description of bound states in a thermal medium has several interesting aspects. When 
the bound state moves with a moderate speed with respect to the medium the resulting EFT is quite similar to the 
one developed for the bound state at rest. We have taken into account the suitable modifications in Section Hill for 
the hydrogen atom in the cases T <^ 1/r and T ^ 1/r. However, when the speed is close to 1, one has to consider two 
well separated scales, and r_, defined in Eq. ([S]), and in the corresponding EFT one has collinear as well as soft 



This is so for generic 9, meaning 9 no ■k/2. For 9 ~ vr/2, we obtain m?j^{v, 9) m^/{l — ti^). 
^ This is so for generic 9, meaning 9 7r/2. For 9 ~ 7r/2, we obtain tic ~ Vl — ae'^ . 



mD{v,9) = mn(0,0) 



h{v,9) 



(91) 



(l~i;2)i/4 ' 




(92) 
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degrees of freedom. The effective temperatures T_|_ and T_ can be in two different energy ranges and in Section HVl we 
have considered two specific cases: the first one corresponds to T_|_ ~ l/r^r_ ^ E and the second one corresponds 
to T-|_ ^ me ^ 1/r ^ ^ E. Note that in this case large logarithms of T^/T^ appear in the calculation. The 
factorized results displayed in Appendix [B] may be useful for a resummation of these large logarithms. It is reassuring 
that our results for moderate velocities are able to reproduce the ones obtained for the w ~ 1 case. For all the cases 
above we observe that the thermal decay width monotonically decreases with the velocity. This means that the faster 
the bound state moves across the thermal bath the more stable it becomes. 

Finally, in Section |V] we have considered the case T ^ 1/r allowing for light fermion pairs in the thermal bath. 
In atomic physics this state could be the muonic hydrogen in a thermal bath of electrons and positrons, while in 
heavy-ion collisions it may represent heavy quarkonia in the quark-gluon plasma at very high temperatures. We have 
determined how the imaginary and real component of the two-body potential are modified for nonvanishing velocities 
of the bound state with respect to the medium. Regarding the real part of the potential we have reproduced known 
results, and extended them to higher speeds. The imaginary part has been calculated for the first time. Its behavior 
is similar to the one determined for the thermal bath at rest for moderate velocities, but it tends to zero at velocities 
close to 1. This implies that Landau damping [3 [H, [13 is not the relevant mechanism for dissociation of bound 
states from a certain critical velocity Vc on, which has been estimated in the previous section. Screening, as originally 
proposed by Matsui and Satz becomes then the relevant mechanism. Our results for the thermal decay width 
disagree with the qualitative estimate of ref. , and with the more quantitative one of ref. [4^ . We believe that 
the main reason for the discrepancy is due to the fact that the velocity dependence of the interaction is not properly 
taken into account in those works. Note that our results for the imaginary part depend crucially on the use of the 
correct non-equilibrium expression in Eq. (j75p . which leads to Eq. (|90p . 

In the present paper we have paved the way for a more detailed study of the propagation of bound states in a 
thermal medium. We have assumed that the medium is a weakly coupled plasma, moving homogeneously and at a 
constant temperature, therefore our study needs a number of refinements to be realistically applied to HQ states in 
heavy-ion collisions. In that case, one should consider the expansion and cooling of the thermal medium, as well as 
possible anisotropics [ill-Hl- In any case, we expect that the qualitative features we observe, namely that the decay 
width decreases with increasing velocity, and hence that Landau damping ceases to be the relevant mechanism for 
dissociation at a certain critical velocity, will remain true. 
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Appendix A: Computation of A and B from section IIIII 

As a starting point one can use Eqs. (16) and (17) from [l3], which for a bound state in a static thermal bath give 



^^K'?) = :e2 — )+^*(^)) • (Al) 



This equation is obtained by a trivial angular integration because the system is symmetric under space rotations. 
When the bound state moves with respect to the thermal bath, the distribution function depends on the effective 
temperature defined in Eq. ([5|) and has a nontrivial dependence of the angle between k and v. Now if we take into 
account the structure that was shown in Eq. (jisp and define x — cos 9, we have that 



and for the real part 



1 



^Aiq) = - / dx ^^^ . (A3) 



21 



Using Eq. ([15]) and the above equations we obtain that the imaginary and real part of the coefficient B are respectively 
given by 



Svr 



^B{q) = / dxx^ log 



27rrVl - f 2 
\q\{l - vx) 



/ j|<7|(l — vx) 



and 



1 

sftB(g) = — / dxx^ 



\q\ 



A 



(A4) 



(A5) 



As a cross-check, in the v — limit we find that the relation B = is fulfilled, and combining this with the identity 

% = + ' (A6) 

one recovers the results reported in fl3 |. 



Appendix B: Computation of the contribution from ultrasoft photons in section IIVI 



In this appendix we compute the matrix elements of Kij defined in Eq. ()56p in the various integration regions 
identified in Section ITVl 



1. The k+,k- ~ T- region 



The quantities a and b defined in Eq. (I58|) can be computed from Kij as follows 

2 _ D-i v'v^ 



b = 



2D -4 
2 

2D- A' 



2D -A 
D v'v^ 



2D -4 V 



2 -"-^J ' 



and for k^,k- ~ T_ we find that 
K,j{q) = ~2i 



i2n) 



d^k S{k+k^ - kl) ( hkj 



1 



2q 



k+ + k- — ie (fc+ + fc_ — iey 



e"^- - 1 

The first term in the square brackets vanishes by symmetry considerations and we find that 

^Kuiq) = 0, 



27r2 V-D- 4 2 2 2 



1 7 1, , /2r_ 

- + - log TT + log 



and 



11 



1 



67r2 V-C'-4 2 ' 2 



log TT + log 



(Bl) 
(B2) 

(B3) 

(B4) 
(B5) 

(B6) 
(B7) 



2. The case with k+,k- ~ g 



In this region we have that 

K^j{q) = ; 



d^k 



(2^) 



D- 



jS{k+k- 



2T- 

IM 



Uk+ + k^)+ie 



(B8) 
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and it is useful to calculate separately the imaginary part of the integrals with the first and the second terms in the 
brackets. The reason is that the computation of the first term in dimensional regularization is technically difficult, 
while the second one is quite straightforward. 

We first compute the imaginary part of the term lineal in T_ with a cut-off to separate the region k+ ^ q from the 
region with fc_ ~ q{T-/T+). Thus we consider a cut-off A such that 9 » A » q{T^/T+) and we obtain 



POO 

QKu{q) = 2T_ / dk+ 
Jo 



lOf 



A 



k+ + 2q + ie 



and 



log 



A 



k+ + 2q — ie 



log 



2q- 



le 



log 



k^—2q—ie 



2q — ie 



(B9) 



V V-' 



Kij{q)=0. 



(BIO) 



Then, the remaining terms are computed in dimensional regularization. Summing all the terms we obtain that 



^Ku{q) = — 



1 



D-4: 

1 



+ + ^- ologTT + logl ^ ) ) - IZ 



+ 2T_ 



f 

Jo 



D 



dkA 



- - -logTT + log 



47r 



lo: 



k+-\-2q—ie 



k-^- + 2q + ie k+ + 2q — ie 



and 



V V-' 



V V-' 



^Kij{q) = 

QK,,{q) = 



2n 

q 

67r2 



67r ' 



1 



loe 



A 



k+ — 2q + ie 



D 



log 



kj^—2q—ie 



k-\- — 2q — ie 



(Ell) 



(B12) 

(B13) 

(B14) 



3. The case with k+ ^ q, k- q{T./T+ 



In this region we have that 

Ki,{q) = 2i 



d^k 5{k+k_ - kl) 



. f d'^k 
7 W)^ 



k+ 



+ 



hk, 
k^ 



+ ie 



(B15) 



For simplicity, we compute the imaginary part using a cut-off (as we did in the previous subsection) and the real part 
using dimensional regularization. We obtain that the real and imaginary parts of the trace of Kij are, respectively, 
given by 



^Ku{q) = - 



1 1 7 1, 

^+2 + 2-2^°^^^ + ^°^ 



27r2 



P 



P 



2q 



2g, 



log 



AT+ 
k+T_ 



+ 2 log 



where P stands for principal value. Moreover, we find that 



V V- 



■^Kij{q) = 0, 



0. 



(B16) 
(B17) 

(B18) 

(B19) 



Appendix C: Interaction with coUinear photons 



In this appendix we give the detailed form of the part of the EFT's Lagrangians that describes the interaction of 
electrons with collinear photons. General remarks about the matching are given in the corresponding subsections. 
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1. pNRQED Lagrangian for T+ ~ 1/r and 1/r > T_ » 



In this case the pNRQED Lagrangian has the form 



SCpNRQED = Cl 



+ C4 



+ C6 



+ C2 



me (jic*) 



a/3 



me ('T-c') (n9) 



+ C3 



{nd) {nd) {nd) {nd) 



me (nd) (nd) ~^ '^^ m. 



indy [nd) 



+ 



ml 

iCB 

ml 



{ndy (nd) 
^n'^n'^F^, n'^n^Fo^p 



+ 



n''n''F^^ n°'n'^ndF, 



(nd) (n9)2 

a/9 



(nd) {ndy 



(nd) (nd) 



n I n ^^ n^'n-'F,,,n'^n^F^p / 
DaV' - — r—TTr'- , „, ^ 



+ 



(nd) 



(nd) (nd) (nd) (nd) 



(nd) 



(nd) (nd) 



+ 



+ — t(A^ 



(nd) (na)2 
(nd) (nd) 



+ ■ 



(nd) _ 
n'^n^'F,,^ n"F, 



(nd) (nd) 



n^F^j n-'djF,, 



+ 



(nd) (na)2 
(nd) (nd) 



(nd) (nd) 

D,tP 



Tl Qj Fjj^j Tl Fjyj 



(n9)2 (nd) 



+ 



(nd) (nd) 



D,i, (Cl) 



The ip field is the field of the electron in pNRQED (we are using the form of pNRQED shown in Eq. (3) of [21|), Di 
is the covariant derivative (containing ultrasoft photons only) that acts on the field of the electron and Latin indices 
stand always for transverse components. 

The power counting for this Lagrangian works as follows. The leading order terms are the ones proportional to ci 
and C2; note that the two terms in n^F^^ are not of the same order of magnitude as ndAi ^ difiA. Taking only the 
leading term of i^^^, in this two terms they have the form 



Cl 



rUe 



C2- 



me 



(C2) 



we consider that all the components of are of the same approximate size (this may not be true in some specific 
gauges). Starting from (jC2p each term that has an additional covariant derivative acting on an electron field is 

and each 



suppressed by an order a, each transverse derivative acting on suppresses this term by an order 



l+v 
3/2 



(nd) acting on a field suppresses the term by an order j:^- If we consider that ^ a ^ V^^J ^^"^ terms 
shown in Eq. (|Cip provide the complete list of operators at the order ra^ct' . In practice only the two first terms 
contribute to our calculation. The Wilson coefficients are fixed by the matching calculation schematically shown 
in (El). We obtain 



Cl = C2 = C4 = Cg = 



2 ' 



C3 



C8 



C7 



C5 



C6 = 
^2 



cio = -y 



(C3) 
(C4) 
(C5) 



2. NRQED Lagrangian for T+ ~ me and l/r > T_ > £ 



The full expression for the part of the NRQED Lagrangian that deals with the interaction of electrons with coUinear 
photons is 



ij^i^n'^F^in-'F,, , ^p'fi, n'^n-'F^^ n^n^ F,,p , i/>t^ 
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(nd) (nd) (nd) (nd) 



n^'n'-F^, n"Fa, , n'^F^, n^n^F^/a 



(nd) (nd) (nd) (nd) 
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+ C8 
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The power counting for this expression is the same as in Eq. (|Cip . However, we are considering now higher velocities, 
and hence the relative size of and a differs from the previous case. For ^ ^ ( j we have listed 

above all the operators up to the order of rriea^. In practice, only the operators proportional to ci, C2, cn and C13 
contribute to our calculation. The Wilson coefficients are fixed by the matching calculation sketched in (|69|. We 
obtain 
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